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Abstract. Our object of study is non smooth vector fields on R^. We 
apply the techniques of geometric singular perturbations in non smooth 
vector fields after regularization and a blow— up. In this way we are able 
to bring out some results that bridge the space between non— smooth 
dynamical systems presenting typical singularities and singularly per- 
turbed smooth systems. 



1. Introduction 

This work fits within the geometric study of Singular Perturbation Pro- 
blems expressed by vector fields on M^. We study the phase portraits of 
certain non— smooth planar vector fields having a curve S as the discon- 
tinuity set. We present some results in the framework developed by So- 
tomayor and Teixeira in [10] (and extended in [9]) and establish a bridge 
between those systems and the fundamental role played by the Geometric 
Singular Perturbation (abbreviated by GSP) Theory. This transition was 
introduced in papers like [2] and [7] , in dimensions 2 and 3 respectively. Re- 
sults in this context can be found in [8]. We deal with non— smooth vector 
fields presenting structurally unstable configurations and we prove that these 
structurally unstable configurations are carried over the GSP Problem asso- 
ciated. Some good surveys about GSP Theory are [3] and 0], among others. 

Let ^ C be an open set and T, C U given by S = /~^(0), where 
/ : l/( — > M is a smooth function having G M as a regular value (i.e. 
V f{p) / 0, for any p G /~^(0)). Clearly S is the separating boundary of 
the regions S+ = {g e U\f{q) > 0} and S_ = {g G U\fiq) < 0}. We 
can assume that S is represented, locally around a point q = {x,y), by the 
function f{x,y) = x. 

Designate by the space of C— vector fields on a compact set K dU 
endowed with the C^— topology with r > 1 large enough or r = oo. Call 
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= Qr{K, f) the space of vector fields Z : K\T, ^R'^ such that 
^(^ „\ - j ^i^^y)^ for ix,y)e^+, 

where X = {fi,gi), Y = (/2,5'2) are in X^'. The trajectories of Z are solu- 
tions of (7 = Z{q), which has, in general, discontinuous right— hand side. 

In what follows we will use the notation 

X.fip) = {Vf{p),X{p)) and Y.fip) = {Vf{p),Y{p)) . 
We distinguish the following regions on the discontinuity set S : 

► Si C S is the sewing region if {X.f){Y.f) > on Si . 

► S2 C S is the escaping region if (X.f) > and (Y.f) < on S2. 

► S3 C S is the sliding region if (X.f) < and (Y.f) > on S3. 

Consider Z G ff^ . The sliding vector field associated to Z is the vector 
field Z^ tangent to S3 and defined at (7 G S3 by Z^{q) = m — q with m being 
the point of the segment joining q + X{q) and q + Y{q) such that m — g' is 
tangent to S3 (see Figured]). It is clear that if g € S3 then (7 G S2 for —Z 
and then we can define the escaping vector field on S2 associated to Z by 
Z^ = —{—Zy. The sewing vector field associated to Z is the vector field 
Z^ defined in g G Si as an arbitrary convex combination of X{q) and Y{q), 
i.e., Z^{q) = \X{q) + (1 — X)Y{q) where A G [0, 1]. In what follows we use 
the notation Z^ for all these cases. 

Let il*" = Q'''{K,f) be the space of vector fields Z : K ^ M."^ such that 
Z G O^' and Z{q) = Z^{q) for all g G S. We write Z = {X,Y), which we 
will accept to be multivalued in points of S. The basic results of differential 
equations, in this context, were stated by Filippov in [5]. Related theories 
can be found in P fTOl [12] . 




Figure 1. Filippov's convention. 

An approximation of the non— smooth vector field Z = (X, Y) by a 
1— parameter family of smooth vector fields is called an e— regularization 
of Z. We give the details about this process in Section [3l A transition func- 
tion is used to average X and Y in order to get a family of smooth vector 
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fields that approximates Z. The main goal of this process is to deduce certain 
dynamical properties of the non— smooth dynamical system (abbreviated by 
NSDS) from the regularized system. The regularization process developed 
by Sotomayor and Teixeira produces a singular problem for which the dis- 
continuous set is a center manifold. Via a blow up we establish a bridge 
between NSDS and the geometric singular perturbation theory. 

Roughly speaking, the main results of this paper are the following: 

Theorem 1. Consider Z{x,y) = Zx{x,y) = {X{x,y),Yx{x,y)) £ Q"^ , where 
A G (—1,1) C M, a non— smooth planar vector field and T, identified with 
the y—axis. Let the trajectories of X be transverse to S and Yq presenting 
either a hyperbolic saddle q £ T, or a hyperbolic focus q £ or a T, — cusp 
point q. Then there exists a singular perturbation problem 

(1) P = a{r, p, y, X) , y' = rP{r, p, y, A) , 

with r > 0, p £ (0,7r), y £ T, and a and (3 of class such that the 
unfolding of ([1]) produces the same topological behaviors as the unfolding of 
the corresponding normal forms of Z\ presented in subsection 15. 

Theorem 2. Consider Z{x,y) = Z^{x,y) = {X^{x,y),Y^{x,y)) £ Q'^ , 
where either p = \ £ R or p = {X,e) £ E?, a non— smooth planar vec- 
tor field and S identified with the y—axis. Consider that q = (xq,yq) £ S 
is a T,—fold point of both X^ and Y^ when p = or p = {0,0). Then there 
exists a singular perturbation problem 

(2) p' = a{r, p, y, X) , y' = rp{r, p, y, A) , 

with r > , p £ (0, tt) , y £ and a and /? of class C such that the following 
statements holds: 

(a) : For all small neighborhood UofqinT, the region (S2US3)n(C/ — 
{yq}) is homeomorphic to the slow manifold a(0, p, y,X) = of ^ 
where y £ {U - {yq}). 

(b) : The vector field Z^ , on (S2 U S3) n (C/ — {yq}), and the reduced 
problem of ([2]), with y £ (U — {yq}), are topologically equivalent. 

(c) : The slow manifold a{0, p,y,0) = of ([2]), where y = yq, has 
just an horizontal component, i.e., a{0, p,yq,0) = can be identified 
with {{p,y)\p £ (0, vr), y = yq}. Moreover, this configuration is 
structurally unstable. 

The unfolding of ([2]) produces the same topological behaviors as the un- 
folding of the corresponding normal forms of Z\ presented in Table[l\and in 

m- 

Observe that Theorem [2] generalize the Theorem 1.1 of p], because here 
we allow that X.f{q) =Y.f{q)=0. 

The paper is organized as follows: in Section [2] we give the basic theory 
about Non— Smooth Vector Fields on the Plane, in Section [3] we give the 
theory about the regularization process, in Section [4] we present the GSP 
Theory, in Section [5] we present the singularities treated in Theorem [T] and 
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give its normal forms, in Section [6] we present the singularities treated in 
Theorem [2] and give its normal forms and in Section [7] we prove Theorems 
Hand El 

2. Preliminaries 

We say that (7 G S is a T,— regular point if 

(i) X.f{q)Y.f{q) > or 

(ii) X.f{q)Y.f{q) < and Z^{q) 7^ (that is g G E2 U E3 and it is not 
a singular point of Z^). 

The points of S which are not S— regular are called T,— singular. We dis- 
tinguish two subsets in the set of E— singular points: E* and S^. Any q £ 
is called a pseudo equilibrium of Z and it is characterized by Z^{q) = 0. Any 
g G S* is called a tangential singularity and is characterized by Z^{q) ^ 
and X.f{q)Y.f{q) = (g is a contact point of Z^). 

A tangential singularity g G S* is a Tj—fold point of X if X.f{q) = but 
^'^■fil) = ^■{^■f){(l) 7^ 0- Moreover, g G S is a visible (resp. invisible) 
^-fold point of X if X.f{q) = and X^.f{q) > (resp. ^^./(g) < 0). 
We say that g G S* is a T,-cusp point of X if X.f{q) = 0, X'^f{q) = but 
X^/(g) 7^ 0. Moreover, g G S is a natural (resp. inverse) T,—cusp point of 
X if X.fiq) = 0, X2./(g) = and X^.fiq) > (resp. X^.f{q) < 0). 

A pseudo equilibrium g G is a T,— saddle provided one of the following 
condition is satisfied: (i) g G S2 and q is an attractor for Z^ or (ii) g G S3 
and g is a repeller for Z^. A pseudo equilibrium g G is a Ti—repeller 
(resp. Y, — attractor) provided q G II2 (resp. q G S3) and q is a repeller 
(resp. attractor) equilibrium point for Z^. 

3. Regularization 

In this section we present the concept of e— regularization of non— smooth 
vector fields. It was introduced by Sotomayor and Teixeira in [10]. The 
regularization gives the mathematical tool to study the stability of these 
systems, according with the program introduced by Peixoto. The method 
consists in the analysis of the regularized vector field which is a smooth 
approximation of the non— smooth vector field. Using this process we get a 
1— parameter family of vector fields G such that for each eg > fixed 
it satisfies that: 

(i) is equal to X in all points of S+ whose distance to S is bigger 
than eo; 

(ii) Z^Q is equal to Y in all points of S_ whose distance to S is bigger 
than eo- 

Definition 1. A C°° function ip : M — > M is a transition function if(p{x) = 
— 1 for X ^ —1, fix) = 1 for x ^ 1 and (p'{x) > if x £ (—1,1). The 
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e— regularization of Z = {X,Y) is the 1— parameter family G X^' given by 

z.(,).(i + ^)x(,)+(i-^)r(,). 

with ifeix) = ip{x/e), for e > 0. 

4. Singular Perturbations 

Definition 2. Let U O be an open subset and take e ^ 0. A singular 
perturbation problem in U (SP— Problem) is a differential system which can 
be written like 

(3) x' = dx/dr = l{x,y, e), y' = dy/dr = em{x,y,e) 
or equivalently, after the time re— scaling t = er 

(4) ex = edx/dt = l{x,y,e), y = dy/dt = m{x,y,e), 

with (x, y) £ U and I, m smooth in all variables. 

The understanding of the phase portrait of the vector field associated 
to a SP— problem is the main goal of the geometric singular perturbation 
theory (GSP— theory). The techniques of GSP— theory can be used to obtain 
information on the dynamics of ^ for small values of e > 0, mainly in 
searching limit cycles. System ^ is called the fast system, and ^ the slow 
system of SP— problem. Observe that for e > the phase portraits of the 
fast and the slow systems coincide. For e = 0, let S be the set 

S = {{x,y):f{x,y,0) = 0} 

of all singular points of ([3]). We call S the slow manifold of the singular 
perturbation problem and it is important to notice that equation Q defines 
a dynamical system, on S, called the reduced problem: 

f{x,y,0) = 0, y = g{x,y,0). 

Combining results on the dynamics of these two limiting problems, with 
e = 0, one obtains information on the dynamics of for small values of e. 
We refer to [4] for an introduction to the general theory of singular pertur- 
bations. Related problems can be seen in [l], [3] and [llj . 



5. Boundary Bifurcations 

Consider Z = {X,Y). In this section we assume that the trajectories 
of the smooth vector field X is transversal to S and that Y has either a 
hyperbolic saddle or a hyperbolic focus or S— cusp point in S. This confi- 
guration is clearly structurally unstable. We present here its normal forms 
and unfoldings. 
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5.1. Codimension One Normal Forms. Take S as the y— axis, i.e., 
f{x,y) = X and consider the parameter A G (—1, 1). 

• Regular— saddle: Assume that X is transversal to S and that Y has 
a hyperbohc saddle in S. The following normal form generically 
unfolds this configuration. 



Z{x,y) = Zx{x,y) 



Xix,y) -- 
yx{x,y) 



-y 
-X — X 



for {x,y)£j:+, 
for (x, y) £ S_. 



Regular— focus: Assume that X is transversal to S and that Y has a 
hyperbolic focus in S. The following normal form generically unfolds 
this configuration. 



Z{x,y) = Zx{x,y) 



X{x,y) : 
y\{x,y) 



X + y + \ 
—X + y — X 



for {x,y)eT,+ , 
, for 



Regular— cusp: Assume that X is transversal to S and that Y has 
a S— cusp point. The following normal form generically unfolds this 
configuration. 



Z{x,y) = Zxix,y) 



X{x,y) : 

Yxix.y) 



1 



for (x, y) G S^ 
, for (x, y) G S_ 



5.2. Regular— saddle Bifurcation. Consider the regular— saddle normal 
form given in the previous subsection. The regularized vector field becomes 

1 — X - X .„sl + x + A 

y = + 

where ip{x/e) is the transition function. Making the polar blow up 



(5) X = r cos 9 and 

where 9 = {2p + 7r)/4, we obtain 



r sm f 



sm( 



2 2 ^ 

1 — r cos 9 — X , ^, 1 + r cos 9 + X 
+ V?(cot0) . 
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Remark 1. We use the new variable 9 = {2p + vr)/4, with p G (0, vr), 
because in the coordinates (r, 9) the transition function ip is constant when 
9 G (0, vr/4) U (37r/4, vr). So, in the next Figures in the text, we can draw the 
slow manifold and slow dynamics with p G (0,7r). 

In the blowing up locus r = the fast dynamics is determined by the 
system 



.sin0(^l_l + ^(cot0)ii^) , y' = 0: 



and the slow dynamics on the slow manifold is determined by the reduced 
system 

— —^ + ip{cot9) — — ^ = 0, y = —— + ip{cot9)——. 



We remark that the slow manifold is implicitly defined by(— l+y)/2 + 
ip{cot 9){—l — y)/2 = and do not depends on the parameter A (see Figure 
[2]). Moreover, y{9) defined in this way is such that 

lim y{9) = +oo and lim y{9) = 0. 




Figure 2. In this figure is pictured the slow manifold to the left 
and the case A < at the right. In both we consider p G (0,7r). 



By other hand, the dynamics on the slow manifold depends on A. In fact, 
if either A > or if A = then y / (see Figure [3]) and if A > so y 
has an unique repeller critical point P (see Figure [2]) given implicitly by the 
equation (^(cot 9) = (-1 + A)/(l + A). 

In Figure [2] and in the next ones, double arrow over a curve means that it 
is a trajectory of the fast dynamical system, and simple arrow means that 
it is a trajectory of the one dimensional slow dynamical system. 

5.3. Regular— focus Bifurcation. Consider the regular— focus normal form 
given in Subsection 15.11 The regularized vector field becomes 
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1+X+x+y ..l-\-x-y 

X = ^ + ^{-e) 2 ' 

l-X-x + y ^,^.1 + X + x-y 

y = 2 + ^(7) 2 ■ 

Similarly to the previous case, considering the polar blow— up given in 
(0), where 9 = {2p + 7r)/4, we get 

a /l + A + y + rcos6' l-X-y-rcos9 
r9 = — smt/l \-(p[cot9) 

1 — A + y — rcos^ , ^, 1 + A — y + r cos ^ 
y = 2 + ^{cot9) . 

Putting r = 0, the fast dynamics is determined by the system 

^' = sm^(^i^ + ^(cot.)^i±^), y' = 0; 

and the slow dynamics on the slow manifold is determined by the reduced 
system 

"'"^"'■+^(cote)^i±A±l = 0, y='-4±y+^icot9)'^^-y 



2 ' 2 ' " 2 ' 2 

In this case, the slow manifold depends of the parameter A. In fact, it is 
given implicitly by (—1 — A — y)/2 + (/?(cot 9){—l + A + y)/2 = (see Figure 
H]). The slow manifold y{9), given in the previous equation, satisfies 

lim y(9) = —00 and lim y(9) = —A. 

4 ^ ^ 4 

We give now the dynamics on the slow manifold. If A < so y 7^ (see 
Figure H]), if A > so y has an unique critical point P, given implicitly as 
the solution of (/?(cot 9) = (—1 + A + y)/(l + A — y), which is an attractor 
(see Figure [5]) and if A = so y / (see Figure [5]) . 



BIFURCATIONS OF NON SMOOTH VECTOR FIELDS ON 9 




Figure 4. In this figure is pictured the slow manifold to the left 
and the case A < at the right. In both we consider p e (0,7r). 




Figure 5. In this figure is pictured the case A = to the left and 
the case A > at the right. In both we consider p G (0,7r). 



5.4. Regular— Cusp Bifurcation. Consider the regular— cusp normal form 
given in Subsection 15. 11 The regularized vector field becomes 



l + A-y2 fx\ l-A + ?/2 



y = LP 

■ € 



2 ^ Ve 

X 



Making the polar blow— up given in ([5]), where 6 = {2p + 7r)/4, we get 



sin( 



((^(cot0)(-l + A-y2)_i_A + y2) 



ij = ip{coi6). 
Putting r = the fast dynamics is determined by the system 

e' = ^((/.(cot + A - y2) _ 1 _ A + ^ y ^ . 
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and the slow dynamics on the slow manifold is determined by the reduced 
system 

ip{cote){-l + A-y^)-l-A + y^ = 0, y = (p{cot0). 

Observe that the slow manifold depends of the parameter A. We can 
obtain the explicit form. In fact, the slow manifold is composed by two 
branches (see Figure [6]): 

/ A(l - v?(cot 6*)) + 1 + ip{cot 6) 



(6) ym 



1 - ip{cot 6) 



The slow manifold satisfies the properties: 
(i) lim y\{e) = ±oo; 



4 



(iz) If A < there exists 9* e (7r/4, 37r/4) (respectively, p* e (0, vr)) 
such that y^{6*) = and the slow manifold is not defined for p G 
{p*,tt). For 9 G (7r/4, 0*) there exist homeomorphisms ^± between 
each branch of the slow manifold and M*. That is, for each z G M* 



there exists 9{z) G (7r/4, 6**) such that y^{9{z)) 



z: 



(Hi) If A > there exist homeomorphisms ^± between each branch of 
the slow manifold and M*. That is, for each z G M* there exists 
9{z) G (7r/4,37r/4) such that yi{9{z)) = z. 

In fact, the item (i) is a straightforward calculus. In order to prove the 
item {a) observe Expression ([6]). Let 9* be such that (/?(cot 9*) = (1 + A)(A — 
1). We have, 

y±io*) = o 

and the radical in ([6]) is negative for 9 G {9*,tt). 
We define the maps: 



C± : R* ^{tt/4.,9*) 

z ^ 9{z) = cot"i (^-1 (^^ITT 

Given z G M* if we put i{z) = 9{z) in ^ we get y±{9{z)) = z. Note that 
^± are homeomorphisms. 

The proof of the item [Hi) is analogous. 

The dynamics in the slow manifold is given by y = (p{cot9). So there 
exists a unique critical point p given implicitly as the solution of (^(cot 9p) = 
0. Note that this critical point is a repeller because (^(cot 0) < for 9 < 9p 
and ip{cot 9) > iov 9 > 9p. See Figure [H 
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6. FoLD-FOLD Bifurcations 

In this section we analyze the dynamics of a NSDS around a point q which 
is S— fold of both X and Y. We say that g is a Fold— Fold singularity of 
Z G Q*". We divide the fold— fold singularities in types according with the 
sign of X'^.f{q) and Y'^.f{q): 

(a) Elliptic case: X'^.f{q) > and Y^.f{q) < 0. See Figure [7] (a) . 

(b) Hyperbolic case: X^.f{q) < 0,Y^.f{q) > 0. See Figured (b). 
(c.l) Parabolic visible case: X^.f{q) > 0,Y'^.f{q) > 0. See Figure El 

(c.l). 

(c.2) Parabolic invisible case: X'^.f{q) < 0,Y'^.f{q) < 0. See FigureE 
(c.2). 

Note that, we can define a first return map ■0^ only in the elliptic case. 
Take E as the y— axis, i.e., f{x,y) = x and consider the parameter A € 
(— 1, 1). The generic normal forms of the hyperbolic and parabolic fold— fold 
singularities are given in Tabled! The normal form of the elliptic fold— fold 
singularity is given in Subsection 16. 4[ 

In the next three subsections we study the dynamics of the hyperbolic 
and parabolic fold— fold singularities via geometric singular perturbations. 
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X 



Y 



(i>) 

Figure 7. Fold-fold singularities. 



X 


Y 






< 




(c. 


1) 



V 

X 

(c.2) 



Hyperbolic 


Parabolic visible 


Parabolic invisible 


= (y-A,-l) 
Y(x,v) = 


Xx{x,y) = {y - \ l) 
Y{x,y) = (-y,-l) 


Xx{x,y) = {y-\-l) 
Y{x,y) = (-y,l) 



Table 1. 



6.1. Hyperbolic Case. Consider the normal form of the hyperbolic fold— fold 
singularity given in Tabled) The regularized vector field is given by 



A fx\f—\ + 2y 



y = -1. 
By the polar blow up we get 

rO = sine Q + (p{cot 9) ^ 

y = -1 
where 6 = {2p + it) /A. 

Putting r = the fast dynamics is determined by the system 

^' = sin0Q + V9(cot0)(^^)) , y' = 0; 

and the slow dynamics on the slow manifold is determined by the reduced 
system 

^ + ^(cot^)(A_^)=0, y = -l. 

In this case we obtain the explicit expression for the slow manifold: 

A(l + yp(cotg)) 
2(^(cot^) ■ 

Observe that, the slow manifold y{9) is not defined for 6q such that ip{cot 6q) = 
0. So, for A / 0, y{9) have two branches and satisfies: 
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(a) lim y{9) = — oo for A < and lim y{9) = +00 for A > 0; 

e — e — 

(6) lim y{6) = +00 for A < and lim y{6) = —00 for A > 0. 

(c) For A = the slow manifold is given implicitly by y(p{cot 6) = 0, 
that is, {{9,y)\9 = Sq} U {(0, y) | y = 0} is the slow manifold. 




A > 



Figure 8. Slow manifold depending of the parameter A. 



The dynamics on the slow manifold is given by y = —1. Therefore, do 
not exist critical points. See Figure [HI 

6.2. Parabolic visible case. Consider the normal form of the parabolic 
visible fold— fold singularity given in Table [TJ The regularized vector field is 

By the polar blow up we get 
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sin^( ^ + ^(cot^) 



ij = Lp{coi6). 

where 6 = {2p + 7r)/4. 

Putting r = the fast dynamics is determined by the system 

^' = sin0Q + ^(cot0)(^^)) , y' = 0; 

and the slow dynamics on the slow manifold is determined by the reduced 
system 



, y = 93(cot 9). 



The analysis is similar to the hyperbolic case. In the present case the 
dynamics on the slow manifold is given by y = (/^(cot 9). So, for A = 0, the 
straight line 9 = 9q is composed by critical points. See Figure [H 




A > 

Figure 9. Bifurcation Diagram of the Parabolic Visible 
Fold— Fold Singularity. 
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6.3. Parabolic invisible case. For this case, we get one different topologi- 
cal type of bifurcation. Consider the normal form of the parabolic invisible 
fold— fold singularity given in Table [H The regularized vector field is 



2 ' Ve/ 2 



y = -1. 
By the polar blow up we get 

re = !l|^(A(^(cot0) + A-2y) , 

y = -1. 

where 9 = {2p + 7r)/4. 

Putting r = the fast dynamics is determined by the system 

Sill 

9' = ^{Xip{cote) + X-2y) , y' = 0; 

and the slow dynamics on the slow manifold is determined by the reduced 
system 

sin^Q^(cot0)+(^^)) =0, y = -l. 
We have the explicit expression for the slow manifold in this case: 

y{9) = ^il + ^{cote)). 

The analysis is similar to the previous cases and the bifurcation diagram 
is expressed in Figure [TDJ 

6.4. Elliptic case. In this case, associated with the non— smooth vector 
fields, there exist the first return map ipzip)- Therefore, we need to analyze 
the structural stability of this one dimensional diffeomorphism. 

Consider Z presenting an elliptic fold— fold singularity, f(x,y) = x and 

^ ^ \ Y{x,y) = {y,-l), for(x,y)GS_. 
The regularized vector field is 

/x\ A 2y — A 



2 2 



y = V^lf 



By the polar blow up we get 
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Figure 10. Bifurcation Diagram of the Parabolic Invisible 
Fold— Fold Singularity. 



. J X^jcote) , (\-2y 
sm p ( ^ h ' 



ij = (/9(cot 9) 

where 9 = {2p + 7r)/4. 

Putting r = the fast dynamics is determined by the system 



!/' = 0; 



and the slow dynamics on the slow manifold is determined by the reduced 
system 



2 V 2 

In this case, for A = 0, we only have sewing region on the non— smooth 
manifold S. The explicit expression for the slow manifold is 
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and there exist one critical point which is attractor if A < and repeher if 
A > 0. See Figure ttH 




Differently of the hyperbolic and parabolic cases, the unfolding (jS]) does 
not give the generic unfolding of a non— smooth vector field presenting a 
elliptic fold— fold singularity. 

Let Zx = {Xx,Y) be as in The expression of its first return map is 



(9) V'za (y) = iY° ixx (y) = y- 2A, 

where 7v>>(p) (respectively 7y(p)) is the first return to S of the trajectory 
of Xx (respectively Y) that passes through p. 

Therefore, we conclude that the critical point of ipz^ is not hyperbolic. In 
order to obtain the generic unfolding of this case we need to unfold the first 
return map. So, the unfolding of the elliptic fold— fold singularity depends 
on two parameters. The first. A, is responsible by the displacement of one 
fold along the y— axis and another one, e, for the unfolding of ipZx- 

Consider the e— perturbation of the smooth vector field Y, given by: 
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Ye{x, y) = gl{x{t),y{t))^ + gl{x{t),y{t))^ 

ox oy 

(10) 

So, the flows of X\ and Yf, are: 

<t>^Xy, (^0, yo) = [xq + (yo - A)t + yo + t), 

4>Y, {xq, yo) = [xo + /d ds, yo + ds 
Let t* G M* and h = 2(A - yo) such that 



(11) / 

Jo 



t* 

£ 



gl{s)ds = Q 



and 0^jO,yo) = (0,-yo + 2A) G S. 

Observe that there exist t* as in Equation ([IT]) , because is a elhptic fold 
singularity for Y . 

We suppose that gf{-),i = 1,2 satisfies: 

(a) Qi are C functions for i = 1,2; 
(6) y,./(0,0)=5f(0) = 0; 

(c) y/./(0,0) =5f(0)^fff(0) + <7i(0)^9f(0) / 0; 

(d) /o**5|(s)d5 = (e-2)y + 0(y2). 

The smooth vector fields Xx,Y^ exhibited in (l8|), ([TO]) , respectively, supply 
the unfolding of ipz- 

^z,,M = <Py, ° ^xS^^y) = (1 - - 2A + 2Ae + 0{y^). 
Therefore, the generic unfolding for the non— smooth vector field Z with 
the origin is an elliptic fold— fold singularity is Zx^e = i^x^Yi.) where 



(12) Zx,e{x,y) 



Xxix,y) = {y - X,l), if(rE,y)ES+, 
ys{x,y) = {gi{x,y),g^{x,y)), if {x,y) E S_ 

and the smooth function gf,i = 1,2 satisfies the conditions (a), (6), (c) and 
(d) given previously. 

7. Conclusion 

We note that: if for any q £ T, we have that Xf(q) 7^ or Yf(q) 7^ 
then, by Theorem 1.1 of [8], there exists a singular perturbation problem 
such that the sliding region is homeomorphic to the slow manifold and the 
sliding vector field is topologically equivalent to the reduced problem. This 
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fact is useful in the next two subsections. 



7.1. Proof of Theorem [T]. In face of Theorem 1.1 of [8], this theorem is 
the subject of section \5\ Moreover, as we give the topological behavior of 
the cases A < 0, A = and A > it is easy to construct the bifurcation 
diagram of [H 

7.2. Proof of Theorem [2l In this theorem we extend Theorem 1.1 of 
[8] considering that can exists a point q such that X.f{q) = Y.f{q) = 0, 
X'^.fiq) ^ and Y^.f{q) ^ 0. In this way, g is a S-fold point of both X 
and Y. 

Consider a NSDS Z\ = {X, Y) where A £ M is a parameter. If with 
the variation of A G {—£,e), the following behaviors are observable then we 
consider that Zx presents a bifurcation, where e > and small. Consider 
A"*" G (0, e) and A~ G (— e,0). The behaviors are: 

(i) : A change of stability on S, i.e., where Zx+ has a sliding region S3 
the non— smooth vector field Zx- has an escaping region E2. 

(ii) : A change of stability on yx, i.e., there are components of S such 
that the induced flow on the slow manifold is such that yx+ > and 

_yx- < 0. 

(ill): A change of stability of the S— singularity, i.e., Zx+ presents a 

S— attractor and Zx- presents a S— repeller. 
(iv): A change of orientation on Si (the sewer region), i.e., Zx+ and 

Zx- presents distinct orientations on Si. 

In face of these previous observations. Theorem [2] follows straightforward 
from section [6l 

Note that, as we give the topological behavior of the cases A < 0, A = 
and A > it is easy to construct the bifurcation diagram of [2] when A G M. 
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